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Abstract In this review we confront the current theoretical understanding of particle ac- 
celeration at relativistic outflows with recent observational results on various source classes 
thought to involve such outflows, e.g. gamma-ray bursts, active galactic nuclei, and pulsar 
wind nebulae. We highlight the possible contributions of these sources to ultra-high-energy 
cosmic rays. 
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1 Introduction 

High-energy astrophysical phenomena stem from the generation of powerful flows emanat- 
ing from supernova explosions, gamma-ray bursts (GRB), from ejections in the environment 
of black holes or neutron stars that lead to the formation of very strong shocks, at which par- 
ticle acceleration takes place. The new developments in these issues, especially for relativis- 
tic shocks, are based on the interdependence between the shock structure, the generation of 
supra-thermal particles and the generation of turbulence. It is thought, and numerical simula- 
tions support that view, that the penetration of supra-thermal particles in the shock precursor 
generates magnetic turbulence which in turn provides the scattering process needed for par- 
ticle acceleration through the Fermi process. This successful development, first elaborated 
for supernova remnants (SNR), inspired similar investigations for the termination shock of 
GRBs. However, in ultra-relativistic shocks, difficulties arise with the transverse magnetic 
field that places a limitation to particle penetration upstream and that drags particles in the 
downstream flow and makes shock recrossing difficult. It turns out that only sufficiently 
fast micro-turbulence can make the Fermi process operative, as demonstrated by recent nu- 
merical simulations. Following a review of the main observational results on GRBs, active 
galactic nuclei (AGN), and pulsar wind nebulae (PWN), these points are briefly discussed 
and astrophysical consequences are drawn. We describe the role relativistic shocks inside 
relativistic flows, e.g. the internal shocks of the prompt-emission stage of GRBs, may play 
in the generation of ultra-high-energy cosmic rays (UHECR). Noting that the energy re- 
quired for supplying sub-GZK UHECR is huge compared with the available energy budget, 
we also discuss other sources, such as AGN and young pulsars, that may contribute to the 
flux at ultra- high energies, all the more so if the composition is enriched in heavy nuclei, as 
suggested by recent experimental results. 


2 Gamma Ray Bursts 

The first GRB was observed by one of the Vela satellites monitoring for the Nuclear Test 
Ban Treaty in 1967, but the (unexpected) astronomical results were not declassified and 
published for another six years (Klebesadel et al. 1973). For many years the nature of GRBs 
was unknown, since the distance scale was completely unknown. Beginning about 20 years 
ago, the cosmological spatial distribution of GRBs was strongly hinted at due to the very 
isotropic distribution on the sky of GRBs localized by CGRO. The wealth of detailed infor- 
mation garnered in the last seven years by Swift (Gehrels et al. 2004) has taken the study 
of GRBs to the next level, and indeed the current situation is in some sense more confusing 
than our naive pr e-Swift picture (e.g. Gehrels et al. 2009). 

In this subsection we (i) review the basic properties of the two main types of GRBs, 
long and short, (ii) look at long GRBs in more detail, (iii) review the brief history of short 
GRBs and the difficulties entailed in their study, (iv) provide an overview of the acceleration 
processes for GRB jets, and (v) conclude with a recent results on high energy emission 
observed by Fermi. 
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Fig. 1 Representative 
broad-band v F v spectra (Gehrels 
et al. 2009) of a 1GRB (910503) 
(Kaneko et al. 2007) and a sGRB 
(980425) (Kaneko et al. 2008) 
along with the Crab pulsar nebula 
(Kuiper et al. 2001) and Cyg X-l 
(McConnell et al. 2002) 



energy [keV] 


2.1 GRB Properties 

GRBs come in two kinds, long and short, where the dividing line between the two is ~2 s 
(Kouveliotou et al. 1993). Long GRBs (lGRBs) are thought to be due to the collapse of a 
massive star, while short GRBs (sGRBs) are inferred to be neutron star - neutron star (NS- 
NS) mergers. A further division can be made spectrally according to their hardness ratio 
(i.e., ratio of high to low energies). The redshift range is from about 0.2 to 2 for sGRBs, 
with a mean of about 0.4. For lGRBs the range is between about 0.009 and 8.2, with a mean 
of about 2.3. The typical energy release is ~10 49 -10 5 ° erg for sGRBs and ~10 5 °-10 51 erg 
for lGRBs. These ranges are based on observed isotropic-equivalent energies of ~10 51 erg 
for sGRBs and ~10 53 erg for lGRBs, and estimates for jet beaming for each class, 0 f ~ 5° 
for lGRBs and 0j ~ 15° for sGRBs (Burrows et al. 2006; Grupe et al. 2006). Beaming angles 
for sGRBs are still highly uncertain. The corresponding beaming factors f h = \ — cos 9, ~ 
0j 2 /2 are roughly 1/300 for lGRBs and 1/30 for sGRBs. The sGRBs have weaker X-ray 
afterglows, a mean value of ~7 x 10 -1 ° ergcm _2 s _1 versus ~3 x 10 -9 erg cm -2 s -1 for 
lGRBs. Figure 1 shows spectra for several representative GRBs, and two other high-energy 
sources, the Crab nebula and Cyg X-l. 

2.2 LGRBs 

The BeppoSAX mission made the critical discovery of X-ray afterglows of long bursts (Costa 
et al. 1997). With the accompanying discoveries by ground-based telescopes of optical (van 
Paradijs et al. 1997) and radio (Frail et al. 1997) afterglows, lGRBs were found to emanate 
from star forming regions in host galaxies at typical distance of z — 1-2. BeppoSAX and 
the following HETE-2 mission also found evidence of associations of GRBs with Type Ic 
SNe. This supported the growing evidence that lGRBs are caused by “collapsars” where the 
central core of a massive star collapses to a black hole (MacFadyen and Woosley 1999). 

LGRBs are incredibly bright. A typical galaxy at a redshift of only z = 3 is fainter than 
m — 27. Multiwavelength observations of the current record holder, GRB 090432 (at z — 8), 


<£) Springer 


312 


A. Bykov et al. 


Table 1 High z GRBs 

z /||j(Gyr) GRB Optical brightness 


8.3 

13.0 

090423 

O 

<N 

II 

@ 

20 min 

6.7 

12.8 

080813 

K = 19 

@ 

10 min 

6.29 

12.8 

050904 

7= 18 

@ 

3 h 

5.6 

12.6 

060927 

7 = 16 

@ 

2 min 

5.3 

12.6 

050814 

K = 18 

@ 

23 h 

5.11 

12.5 

060522 

7? = 21 

@ 

1.5 h 


are providing information about the universe at a time when it was only about 4 % of its 
current age, and shed light on the process of reionization in the early universe (Tanvir et al. 
2009; Salvaterra et al. 2009). The highest redshift GRBs are seen to have high luminosity, 
resulting in fluxes well above the detection threshold. Such bursts are also strong at other 
wavelengths. Table 1 presents optical data for the highest redshift GRBs observed to date, 
where the look-back time t LB (Gyr) is given in column 2. 

2.3 SGRBs 

At Swiff s launch, the greatest mystery of GRB astronomy was the nature of short-duration, 
hard-spectrum bursts. Although more than 50 lGRBs had afterglow detections, no afterglow 
had been found for any sGRB. Swift provided the first sGRB X-ray afterglow localization 
with GRB 050509B and HETE-2 enabled the first optical afterglow detection with GRB 
050709. These two bursts, plus Swift observations of GRB 050824, led to a breakthrough 
in our understanding (Gehrels et al. 2005; Bloom et al. 2006; Fox et al. 2005; Villasenor et 
al. 2005; Hjorth et al. 2005; Barthelmy et al. 2005; Berger et al. 2005) of sGRBs. BAT has 
now detected 60 sGRBs, most of which have XRT detections, and about one third of which 
have host identifications or redshifts (an additional two have been detected by HETE-2, one 
by INTEGRAL, and two by Fermi/LAT). We now have ~50 sGRB localizations. 

In stark contrast to long bursts, the evidence to date on short bursts is that they can 
originate from regions with low star formation rate. GRB 050509B and 050724 were from 
elliptical galaxies with low current star formation rates while GRB 050709 was from a 
region of a star forming galaxy with no nebulosity or evidence of recent star formation ac- 
tivity in that location. Recent HST observations of locations of sGRBs in their hosts reveal 
that sGRBs trace the light distribution of their hosts while lGRBs are concentrated in the 
brightest regions (Fong et al. 2010). SGRBs are also different from lGRBs in that accom- 
panying supernovae are not detected for nearby events (Bloom et al. 2006; Fox et al. 2005; 
Hjorth et al. 2005). Taken together, these results support the interpretation that short bursts 
are associated with an old stellar population, and may arise from mergers of compact bina- 
ries [i.e., double neutron star or neutron star-black hole (NS-BH) binaries]. 

2.4 GeV Emission 

Fermi was launched into low-Earth orbit in June 2008 and has two primary high-energy de- 
tectors: the Large Area Telescope (LAT) which operates between 20 MeV and > 300 GeV, 
and the Gamma- Ray Burst Monitor (GBM) which operates between 8 keV and 40 MeV. So 
far the LAT has detected 24 GRBs; two were sGRBs, and nine showed extended emission. 
The emission from GRB 090902B included a 34 GeV photon. One of the most lu mi nous 
to date has been GRB 0809 16C (Abdo et al. 2009a) at a redshift of 4.35. It had extended 
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emission (18 min) and exhibited a lag in LAT energies with respect to GBM. GRB 090510 
is unique in being the only short burst with a known redshift (0.903) showing GeV emission. 
The lack of detectable time delay between specific peaks in the light curves of GRB 0905 10 
at different energies leads to strong constraints on Lorentz invariance (Abdo et al. 2009b). 
Recent theoretical work (Kumar and Bamiol Duran 2010) on the Fermi/LAT detected GRBs 
suggests that these may represent unusually powerful explosions with Lorentz factors >10 3 
in which the entire progenitor is obliterated. The simplest model, namely an external shock 
with synchrotron emission, can be used to take the early values (at ~10-10 2 s) of the ob- 
served high-energy emission and successfully predict the much later values of the optical 
and X-ray afterglow (at ~10 5 -10 6 s). 

2.5 GRB Summary 

Recent progress in GRB research has been strongly motivated by observational discoveries. 
To date, Swift has detected about 600 GRBs, Fermi /LAT 24. High redshift GRBs are illumi- 
nating the properties of the high- ^ universe and probing into the era of re-ionization. Swift 
finds sGRBs in different environments than lGRBs; also sGRBs are not accompanied by su- 
pernovae. The accumulating evidence provides support for the NS-NS merger model. Many 
GRBs have delayed onset of GeV emission, and more have extended high energy emission. 
Interesting constraints on the Lorentz factors associated with outflow, and Lorentz invari- 
ance violation, also come from the synergism between GeV and lower energy observations. 
It is not currently know for certain whether GRB jets are made primarily of baryons or 
Poynting flux, but momentum for the latter idea appears to be gaining strength. 

2.6 Jet Launching Processes 

There are currently two primary lines of thought regarding the creation and propagation of 
jets in GRBs. 

The baryonic jet model, whose roots can be traced back to the elegant analytical solutions 
of a relativistic blast wave by Blandford and McKee (1976), posits that a jet containing about 
a Jupiter’s mass worth of gas, - — 0 .00 1 M 0 , is somehow launched near the BH created by the 
collapsar with a Lorentz factor r — 10-20 (Zhang et al. 2003). The jet propagates through 
the dense stellar envelope of the progenitor star where it is focussed and compressionally 
heated. After breaking free of the stellar surface, the thermal energy of the compressed 
jet is translated into bulk kinetic motion, further accelerating the jet to a Lorentz factor of 
several hundred. Subsequent deceleration by the circumstellar medium of the jet, which is 
then idealized as being the fragment of a relativistic shell so that the Blandford & McKee 
formalism can be brought to bear (Sari et al. 1998), can then be used to infer jet beaming 
angles from putative achromatic “breaks” in the GRB decay light curves (Frail et al. 2001). 

The Poynting flux jet model, which has been gaining momentum in recent years, traces 
its roots back to Blandford and Znajek (1977) who considered the electromagnetic (EM) 
extraction of energy from within the ergosphere of a Kerr BH. In the last few years workers 
have developed sophisticated numerical codes that calculate the 3D evolution of gas and EM 
fields from the inner edge of accretion disks onto spinning BHs, taking into account both 
general relativity and magnetohydrodynamics (McKinney 2005; McKinney and Narayan 
2007a, 2007b). The accretion disk inner edge provides a natural collimating surface. These 
workers find a baryonic zone-of-exclusion within the jet, which effectively suppresses any 
baryonic component. For numerical stability the numerical codes need to have some mass 
within each grid point, therefore a small trace amount of matter is constantly added within 
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grid points where the EM fields try to exclude it. These are called “floor” models. Previous 
studies that attempted to place constraints on the jet Lorentz factors, which were based on 
the baryonic jet assumption, have been called into question. The ramifications on the jet 
dynamics of having a predominantly Poynting flux jet have not been developed yet in any 
detail. 


3 Nonthermal Particle Acceleration in GRBs: Challenges and Perspectives 

Physical mechanisms of prompt emission in GRBs are still to be established. There are fun- 
damental questions of the powerful relativistic outflow composition and matter vs. magnetic 
field dominance to be addressed with both observations and advanced models. The relativis- 
tic outflows may be different for two types of GRB progenitors under consideration. Long 
and soft GRBs are most likely connected to supemovae, while the short and hard GRBs are 
possibly related to compact binary mergers. 

Magnetized outflows from GRB engines of different nature have been studied for 
some decades (see, e.g. Usov 1992; Thompson 1994; Meszaros and Rees 1997; Lyu- 
tikov and Blandford 2004; Zhang and Yan 2011). The outflows converting the rotation 
power of a compact collapsar into a broad band of radiation are likely models of GRBs 
and AGNs. The principal question here is our understanding of the microscopic mech- 
anisms of the conversion of magnetic field energy into non-thermal particles and the 
observed emission — with or without shock formation. Anisotropic striped wind with al- 
ternating magnetic polarity is considered a favorable configuration to convert the mag- 
netic energy into the observed radiation of pulsar wind nebulae (e.g. Kirk et al. 2009; 
Sironi and Spitkovsky 2011). An alternating magnetic field configuration may occur in 
relativistic jets of GRBs. Magnetic field reconnection demonstrated to produce electric 
fields that accelerate particles in the Earth magnetosphere, solar flares, and some laboratory 
plasma configurations (see e.g. Yamada et al. 2010) was also proposed as a plausible GRB 
model by (e.g. Spruit et al. 2001; Lyutikov and Blandford 2004; Giannios and Spruit 2005; 
Zhang and Yan 2011; McKinney and Uzdensky 2012). 

The most elaborated model of the origin of the GRB prompt emission is, by now, the 
relativistic dissipative fireball model. There are a number of alternative particle acceleration 
and radiation processes within the relativistic dissipative fireball paradigm (see e.g. Rees 
and Meszaros 1994; Piran 2004; Meszaros 2006; Mimica and Aloy 2012). Synchrotron and 
inverse Compton radiation in the optically thin regions of the relativistic fireball can be as- 
sociated with the non-thermal electron/positron accelerated either in the internal dissipation 
processes due to shocks or to the flow magnetic field reconnections. Another potential com- 
ponent is the photospheric emission (thermal or non-thermal) that is rather rarely identified 
in the observed GRB emission. The prompt emission light curves and spectra are generally 
in agreement with the internal shock models. However, some potential problems of the sce- 
nario are its efficiency and the lack of a bright photospheric component observed in a few 
GRBs. We shall discuss now some general features of the internal dissipation models with 
emphasis on particle acceleration processes. 

3.1 Energetic Efficiency of Internal Dissipation Models: Shell Collisions in Jets 

An important issue of the internal dissipation scenario where the energy carried out by mul- 
tiple colliding shells of different magnetization is the efficiency of the conversion of the out- 
flow power into the observed radiation. The variability of the central engine (of a timescale 
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Aar) driving the relativistic outflows of the mean Lorentz factor A, either matter or Poynt- 
ing dominated, can be modelled as a collision of energy-containing shells. The model is 
considered to explain the main features of the GRB prompt emission (see e.g. Piran 2004; 
Meszaros 2006), as a vital alternative to the photospheric models of GRBs. The dissipation 
region typically exists at the radii about r<ii SS ~ ct vai r 2 . In the case of matter dominated 
jets (of low magnetization) the inner dissipation occurs in the internal shocks while in the 
electromagnetically (Poynting flux) dominated jets the magnetic field reconnection effects 
are most likely crucial though shocks may also occur. The microphysics of the dissipation 
in relativistic shocks as well as modeling of the magnetic field reconnections are under in- 
tense studies (Yamada et al. 2010; Bykov and Treumann 2011; Sironi and Spitkovsky 2011; 
McKinney and Uzdensky 2012). Realistic models of a jet that would simulate the global 
RMHD dynamics and simultaneously resolve the dissipative microphysical plasma pro- 
cesses at much smaller spatial scales are not feasible at the moment. However, simple mul- 
tiple shell models of the internal dissipation that just parameterize the magnetic field recon- 
nection effects are still rather useful and it is instructive to discuss some of these models. 
A similar approach can be applied to other relativistic outflows, like those of the AGN jets 
(see for a discussion Sect. 4) and of the pulsar wind nebulae. 

To illustrate the effect of the outflow magnetization a = B 2 /An T pc 2 on the en- 
ergy conversion into the observed radiation it is instructive to use a simple two shell 
model (see e.g. Panaitescu et al. 1999; Kumar 1999; Daigne and Mochkovitch 1998; 
Zhang and Yan 2011). The shell collision may result in dissipation of the magnetic en- 
ergy due to reconnection and turbulence cascade. Consider two shells of masses, Lorentz 
factors and magnetization parameters [((mi, A, crj) and (m 2 , A), cr 2 )], respectively, collid- 
ing inelastically with the formation of a merged shell of (mf, A, cr f ). The internal energy 8' 
released in the rest frame of the merged shell (of the Lorentz factor A) is assumed to be 
in the form of either thermal or non-thermal accelerated particles and radiation with nearly 
isotropic distribution in the rest frame. Then, in the observer frame the released energy is 
A E = r { 8 This leads to conversion of some amount of magnetic energy into internal en- 
ergy of the fluid, and then to radiation. Let us envisage a picture where the two shells merge 
with a lower magnetization parameter crf by the end of such an inelastic collision. Energy 
conservation and momentum conservation can be presented as 


A A m 1 + A Am 2 = AA(mi + mi + S') 



A (A) = 1 + 


ir 2 - 1 

2 a 2 


°i 


( 1 ) 


and 

AA m i(l + (J\) + AAmi(l + <y 2 ) = Af/Jf(/ni + m 2 + Y)(l + Of). (2) 

Since the merged shell can not be considered as a cold one anymore 


(P' + p'fW 


( 3 ) 


where p’ { is the proper kinetic energy density (i.e. the internal energy density with the rest 
energy density subtracted), V' is the shell volume, and P' is the pressure in the rest frame of 
the merged shell. Note that to calculate the energy released in the rest frame of the merged 
shell 8' one should keep all of the terms of the order of F~ 2 . The expression for A (A) is 
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Fig. 2 Simulated energy dissipation efficiency defined by Eq. (5) as a function of the shell magnetization 
parameter a for the inelastic collision of two initially cold magnetized shells of equal masses, but of different 
Lorentz factors (T\ = 500 and A = 1000). The key parameter is the final magnetization of the hot merged 
shells Cf was chosen to be 0.1 ( dot-dashed curves), 1.0 ( solid lines) and 10.0 ( dashed lines). The two curves 
of the same line style in the Figure differ by the assumed initial magnetization parameters. The top curve for 
each line style corresponds to the case of equal initial magnetization oq = oq. while the lower curves of each 
type correspond to the fixed aq =0.1. The adiabatic index y = 4/3 was fixed for the matter in the hot merged 
shell 


exact for the case of the transverse magnetic held in the rest frame of a shell, and it accounts 
for the energy of the induced electric fields in the observer frame. 

In general, there are a few distinct components that contribute to the pressure and the 
proper energy density (thermal and nonthermal baryons, leptons, and photons). If for sim- 
plicity we describe these as a single fluid with adiabatic index y , then 


, , yP' 

P> + p=^—, (4) 

Y - 1 

and one can resolve the energy-momentum conservation equations, assuming the simple 
equation of state Eq. (4) to calculate the Lorentz factor of the merged shell r f and the energy 
dissipation efficiency rj of the inelastic collision of the two cold magnetized shells 


r f <5' 

r I P i m I + r iPiini ‘ 


( 5 ) 


In Fig. 2 we show the simulated energy dissipation efficiency i) as a function of the shell 
magnetization parameter a for the inelastic collision of two initially cold magnetized shells 
of equal masses, but of different Lorentz factors. We consider the GRB jet as a generic case 
and therefore choose fj = 500 and A = 1000. The energy dissipation efficiency rj, the mean 
enthalpy of the hot matter in the merged shell, that characterizes the mean Lorentz factor of 
the randomized particles in the rest frame of the hot merged shell (shown in Fig. 3), and the 
Lorentz factor of the hot merged shell in the observer frame, shown in Fig. 4, are derived 
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Fig. 3 Simulated dimensionless 
mean enthalpy that is 
characterizing the mean Lorentz 
factor of the randomized particles 
in the rest frame of a hot merged 
shell. The curves are simulated 
for the same parameter sets as it 
is indicated in Fig. 2 


6 

+ 


CD 


£ 

CD 

£ 



Fig. 4 Simulated Lorentz factor 
7f of a hot merged shell in the 
observer frame. The curves are 
simulated for the same parameter 
sets as it is indicated in Fig. 2 



magnetization 


from the conservation laws and the equation of state. The adiabatic index is fixed to y = 4/3, 
though in more accurate numerical simulations it depends on the Lorentz factor /f. The 
key parameter here is the final magnetization of the hot merged shell ctf that is determined 
by the currently poorly known rate of magnetic field reconnection in the merging shells. 
The case of fast field reconnection (and, therefore, efficient magnetic field dissipation) is 
illustrated by the dot-dashed curves in Fig. 2 where ay = 0.1. The cases of lower magnetic 
field dissipation are illustrated by the final magnetization parameters a t = 1.0 (solid lines) 
and erf = 10.0 (dashed lines). 
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The dependence of the merged shell parameters on the initial magnetization of the collid- 
ing shells is presented by two curves of the same style that differ by the initial magnetization 
parameters. The top curve for each type of line corresponds to colliding shells of equal ini- 
tial magnetization a\ = 02 . The lower curves of each type correspond to the case when the 
fast shell of the Lorentz factor A = 1000 has low initial magnetization os = 0.1. 

It is clearly seen in Fig. 2 that the dissipation efficiency is higher in the case of the ini- 
tially highly magnetized shells with fast magnetic field reconnection resulted in low Of of 
the merged shell (dot-dashed curves). Again, the microscopical model of the reconnection 
rate in the complex flow is still to be done to estimate the crucial value of the final mag- 
netization parameter cr t - of the merged shell (see e.g. Yamada et al. 2010). To investigate 
the problem of shock formation in the internal dissipation scenario, a microscopic model- 
ing of the collisionless shock formation and its structure in a highly magnetized relativistic 
outflow is needed, and that is a truly challenging task (see e.g. Bykov and Treumann 2011; 
Sironi and Spitkovsky 2011). The problem of Fermi acceleration in transverse relativistic 
shocks of different magnetization that are important to describe the external shocks in the 
jets of GRBs and AGNs will be addressed in §6, here we concentrate on the internal dissi- 
pation models. 

As it is seen in Fig. 4, the high Lorentz factors A > 1000 of a hot merged shell in the 
observer frame can be achieved even for the incomplete magnetic field dissipation in the 
merged shell (solid and dashed curves) and, therefore, it can further catch up other slowly 
moving shells and merge with them providing a chain dissipation process. 

The mean Lorentz factors in Fig. 3 derived from the conservation laws do not preclude a 
presence of a non-thermal (e.g., a piece-wise power-law) particle distribution, where some 
minor fraction of particles can reach energies that are by some orders of magnitude larger 
then the derived “thermodynamic” mean Lorentz factor. Now we turn to discuss in brief 
possible particle acceleration processes. 

3.2 Particle Acceleration in the Internal Dissipation Models 

Relativistic turbulence produced by the internal dissipation in shocks and magnetic re- 
connection in relativistic jets result in acceleration processes occurring on both gyro time 
scale and on longer comoving hydrodynamical time scales of the order of / /c. The elec- 
tric fields induced by turbulent motions of plasmas carrying magnetic fields on different 
scales lead to statistical energy gains of the superthermal charged particles and their wide- 
band radiation (see Bykov and Meszaros 1996; Mizuno et al. 2011; Nishikawa et al. 2010; 
Sironi and Spitkovsky 2011; Zhang and Yan 2011; Bykov et al. 2012; Murase et al. 2012). 

For nonrelativistic MHD turbulence the particle energy gain over a turbulent corre- 
lation length (or correlation time) is small, because the induced electric field is smaller 
then the entrained magnetic field. However, the distinctive feature of statistical accel- 
eration in the relativistic MHD turbulence and shocks on larger scales expected in the 
flow-colliding regions, is the possibility of a substantial particle energy gain over one 
correlation scale, because the induced electric fields are no longer small. In this case 
the standard Fokker-Planck approach cannot be used. Instead, it is possible to calculate 
the energy spectra of nonthermal particles within a special integro-differential equation 
which is a generalization of the Fokker-Planck approach (see Bykov and Toptygin 1993; 
Bykov and Meszaros 1996). 

Charged particles interact with a wide spectrum of RMHD fields and an internal shock 
ensemble produced by the colliding shells. In the comoving frame, it is assumed that the fluc- 
tuations on all scales up to ~ A (including the internal shock ensemble) are nearly isotropic 
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(in the latter case, it is enough if they are forward-backward symmetric). The small mean 
free path k of the superthermal particles leads to their isotropy in the frame of the local bulk 
velocity fluctuations. The assumed statistical isotropy of the bulk velocity fluctuations in the 
comoving frame of the wind results then in a nearly isotropical particle distribution, after 
averaging over the ensemble of strong fluctuations on scales ~Z. 

To calculate the spectrum of nonthermal leptons accelerated by an ensemble of inter- 
nal shocks and large-scale plasma motions in the flow-colliding region (averaged over the 
statistical ensemble of large-scale field fluctuations) we use a kinetic equation for the nearly- 
isotropic distribution function N = y 2 F, which takes into account the non-Fokker- Planck 
behavior of the system (see Bykov and Toptygin 1993; Bykov and Meszaros 1996), 

3F(r,$,f) f 00 

' = 6,~(g)+ / d&ZMS-SOAFfcSi.f) 

Ot J - oo 

/a 2 aw 00 

Here £ = 1 n ( y / y, ) , y, is the Lorentz factor of the injected particles, £?,•(§) oc t;cnl 2 is the rate 
of nonthermal particle injection, n is the lepton number density in the local flow comoving 
frame. The kernels of the integral equation (6) determining the spatial and momentum dif- 
fusion are expressed through correlation functions which describe the statistical properties 
of the large scale MHD turbulence and the shock ensemble. Following the renormalization 
approach, the Fourier transforms of the kernels (s) and ( 5 ) are solutions of a tran- 
scendental algebraic system of equations of the form D[\ = d> j2 (Df. Zkf , s). Here s is a 
variable which is the Fourier conjugate of £. Equation (6) and the renormalization equations 
are valid only for particles with sufficiently small mean free paths A(y) « A. 

It is important to note that the solution of Eq. (6) has a universal behavior, only weakly 
dependent on the complicated details of the turbulent system. The stationary solution to 
Eq. (6) with a monoenergetic injection rate Q t has an asymptotical behavior of a power-law 
form, N oc 2, y _a , where a = —0.5 + [2.25 + (O)/)^ -1 (O)] 0 5 , and thus one may take 6 ~ 

(//A) 2 . For conditions typical of a developed RMHD turbulence, the ratio of the rate of the 
scatterings to the acceleration rate is D 1 (0) D 1 (0) < 1 , and for 9 < 1 one obtains a ~ 1 . This 
hard y “ 1 spectral behavior arises because the acceleration time z a ~ l/c ~ a A /c is much 
shorter than the escape time at the relevant energies, r eS c ~ A 2 /k ~ A 2 /(/c) ~ A /( ac ). The 
power needed to produce such a spectrum of nonthermal particles increases oc y max , so it is 
important to understand its temporal evolution. 

In the test particle limit, where the backreaction of the accelerated leptons onto the 
energy-containing bulk motions is small, we have N(y, t ) oc £ ny~ l for y < y*(f), where 
Y*(t) = Yi expa/rj') and 

r* oc//c~a(A/c), (7) 

is the typical hydrodynamical acceleration timescale (see, e.g. Bykov and Toptygin 1993), 
with y, ~ few, a < 1, and the comoving width of the region energized by shocks equal to A. 
From the energy balance equation, when the value y*(t) ~ yim p /m e e^~ l is reached, the 
growth must saturate, and the resulting spectrum consists of two branches. One of these is 
the hard spectrum N(y) ~ £«y _1 > f° r Y < Y*, where 

Yi m p /m e s £ _1 . (8) 

For the typical scales of the considered problem y, ~ 1 and £ ~ 10 -3 , so y* ~ 10 5 (but 
it could be even larger since s ~ 1 for large scale plasma motions). In the regime of 
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very efficient particle acceleration where the backreaction of the accelerated leptons on the 
energy-containing bulk motions is substantial, one should use a non-linear approach. In 
that case the particle distribution experience fast temporal evolution (see e.g. Bykov 2001; 
Ferrand and Marcowith 2010). 

The multiwavelength observations discussed above can be used to constrain the charac- 
teristics of relativistic turbulence and particle spectra evolution. Measurements of polariza- 
tion of prompt gamma-ray emission of GRBs would be very informative. Polarization mea- 
surements provide an additional deep insight into the nature of acceleration/radiation pro- 
cesses. By now only a few instruments are available for this kind of measurement at gamma- 
ray energies. Gamma-ray bursts (GRBs) are very promising candidates for polarimetry due 
to their large flux over limited time intervals, maximizing the available signal-to-noise ratio. 
To date, a few polarization measurements have been reported, all claiming a high degree 
of polarization in the prompt emission of GRBs, but with rather low statistical evidence 
yet. Gotz et al. (2009) used the IBIS telescope aboard the INTEGRAL space laboratory to 
measure the polarization of the prompt gamma-ray emission of the long and bright GRB 
0412 19 A in the 200-800 keV energy band. They found a variable degree of polarization 
ranging from less than 4 % over the first peak to 43 %-25 % for the whole second peak. 
Time-resolved analysis of both peaks indicated a high degree of polarization, and the small 
average polarization in the first peak can be explained by the rapid variations observed in the 
polarization angle and degree. The results by Gotz et al. (2009) are consistent with different 
models for the prompt emission of GRBs at these energies, but they favor synchrotron radi- 
ation from a relativistic outflow with a magnetic field, which is coherent on an angular size 
comparable with the angular size of the emitting region. Recently, Yonetoku et al. (2011) re- 
ported the polarization measurement of the prompt gamma-ray emission of GRB 100826A 
with the Gamma-Ray Burst Polarimeter (GAP) aboard the small solar-sail demonstrator 
IKAROS. They detected the change of polarization angle during the prompt emission, and 
the average polarization degree of 27 ± 1 1 %. Polarization measurements are a powerful tool 
to constrain the GRB radiation mechanisms and the magnetic field structure in the photon 
emitting regions. 


4 Active Galactic Nuclei 

4.1 Studies of Particle Acceleration in AGN Jets 

The relativistic outflows of Active Galactic Nuclei (AGNs) are well-studied particle acceler- 
ators. The AGN jet phenomenon spans many orders of magnitudes: jets originate on sub-pc 
scales close to supermassive black holes with typical masses between a few million and a 
few billion solar masses and can propagate over many hundred kpc to feed giant hot-spot 
complexes and radio lobes. One of the most remarkable properties of jets is that they dis- 
sipate little energy while transporting vast amounts of energy and momentum over large 
distances. However, jets are not dissipation-free neither one sub-pc scales nor on kpc-scales: 
they do accelerate electrons and maybe also protons to high energies. The high-energy par- 
ticles interact with magnetic fields, photon fields, and/or ambient matter and make AGNs 
some of the brightest extragalactic sources of continuum emission across the electromag- 
netic spectmm. There are many reasons to study particle acceleration in AGNs: we would 
like to explain the observed electromagnetic radiation in order to constrain the composi- 
tion and structure of jets. The composition and structure constrain the processes of AGN 
accretion and jet formation, acceleration and collimation. The same studies can enhance 
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our understanding of AGN feedback, i.e. how AGNs interact with their hosts and decisively 
impact their evolution. Studies of AGNs also allow us to study the particle acceleration 
mechanisms. A good (“microscopic”) understanding of the dissipation processes in jets is 
instrumental for addressing the larger questions mentioned before. 

In the following we discuss two particle acceleration sites: the “blazar zone” less than a 
pc away from the supermassive black hole, and jet particle acceleration by kpc-scale jets. 
The processes at work in these two well- separated regions have recently received a lot of 
attention mainly because of spectacular jet images from the VLBA, HST, and Chandra, and 
spectacular multiwavelength observations with the RXTE, Suzaku, and Swift, X-ray, and 
Fermi, H.E.S.S., MAGIC, and VERITAS y-ray telescopes. 

4.2 Studies of Particle Acceleration Processes in the Blazar Zone 

The cores of the “blazar class” of AGNs are bright sources of spatially unresolved continuum 
emission. In the case of BL Lac objects, emission and absorption lines are absent or weak. In 
the case of their more powerful siblings, Flat Spectrum Radio Quasars (FSRQs), the spectra 
exhibit emission and absorption lines, which afford additional diagnostics concerning the 
mass of the central black hole, the accretion power, and the redshift of the source. Blazars 
are sources with jets closely aligned with the line of sight. The relativistic propagation of 
the jet plasma and the close alignment of the jet with the line of sight lead to relativistic 
boosting of the emission. Blazar Spectral Energy Distributions (SEDs) show evidence for 
two broad peaks presumably coming from synchrotron and inverse Compton emission from 
the same electron population. Observations of the BL Lac Mrk 501 revealed emission up 
to 16 TeV (Aharonian et al. 1999), giving direct proof that AGN jets accelerate particles to 
5>> 1 TeV energies. AGNs may accelerate particles to much higher energies, maybe even to 
ultra-high energies (e.g. KachelrieBet al. 2010), see however Lemoine and Waxman (2009) 
for a detailed discussion of acceleration power in these sources. 

Two scenarios are commonly invoked to explain the observed emission: (i) the jet is 
initially Poynting flux dominated and accelerates particles in magnetic reconnection events 
(e.g. Giannios et al. 2010, and references therein); (ii) the jet is either particle energy dom- 
inated right from the start, or, it is initially Poynting flux dominated and some unknown 
mechanism converts the energy from Poynting flux into bulk motion energy, and the particle 
dominated jets subsequently accelerate particles at shocks. 

Blazar studies of particle acceleration benefit from the fact that the broad-band SEDs 
oft he sources evolve on short time scales: X-ray and gamma-ray flares with 1 min dura- 
tions have been reported. It thus is possible to track the temporal evolution of the particle 
energy spectra. Modeling of sequences of broad-band SEDs has shown that the jets are 
highly relativistic with bulk Lorentz factors of ~50 or even higher (e.g. Krawczynski et 
al. 2001). The simplest types of models synchrotron self-Compton models — indicate that 
the particle energy density dominates strongly over the magnetic field energy density in the 
blazar zone (e.g. Krawczynski et al. 2002; Acciari et al. 2011). These results clearly favor 
shock acceleration over magnetic reconnection. The data and theoretical studies (particle- 
in-cell (PIC) simulations and instability analyses) favor mildly relativistic shocks internal to 
the jets over highly relativistic external shocks: the former have faster downstream plasmas 
commensurate with the high inferred bulk Lorentz factors of the emitting plasma, and are 
less susceptible to the strong suppression of the formation of upstream scattering centers by 
even a weak plasma magnetization than their highly relativistic counterparts (e.g. Lemoine 
and Pelletier 2010; Sironi and Spitkovsky 2009). 

X-ray and very high-energy y-ray observatories have recently succeeded to sample the 
temporal evolution of the X-ray and X-ray energy spectra with delicate accuracy. Although 
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observers organized a large number of observation campaigns with the objective to find 
the flux vs. spectral index evolution patterns predicted by simple acceleration theory (e.g. 
Kirk et al. 1998), the observations revealed rather erratic evolutions (e.g. Takahashi et al. 
2000; Garson et al. 2010, and references therein). Unfortunately, observations have not yet 
succeeded to determine unambiguously where exactly the blazar emission originates. For 
the radio galaxy M87 possibly a misaligned blazar the observation of temporally coincident 
radio, X-ray, and y-ray flare indicates that the emission comes from <100 Schwarzschild 
radii of the supermassive black hole (Acciari et al. 2009). 

4.3 Studies of Particle Acceleration Processes of kpc-scale Radio, Optical, and X-Ray Jets 

The VLBA, HST, and Chandra telescopes are delivering images of the kpc-jets of an ever 
increasing number of radio galaxies. The images in the radio, optical, and X-ray bands can 
be used to infer complementary information about the emitting particle populations. Radio 
and optical polarimetry give additional clues about the orientation of the magnetic field in 
the bright jet regions. 

A recent somewhat surprising result was the detection of a large number of X-ray bright 
kpc-jets with Chandra (see the review by Flarris and Krawczynski 2006, and references 
therein). In the case of powerful Fanaroff-Riley (FR) type II sources (like PKS 0637-752) 
the combined radio, optical and X-ray energy spectra showed unambiguous evidence for 
two distinct emission components. Presently two models are commonly invoked to explain 
the “second” X-ray bright component: (i) inverse Compton scattering of Cosmic Microwave 
Background photons (Tavecchio et al. 2000; Celotti et al. 2001), (ii) a second synchrotron 
component (Harris and Krawczynski 2002). The first model requires jet plasma moving with 
large bulk Lorentz factors (C ~ 10) at kpc-distances from the central engine. The relativistic 
motion boosts the CMB photon energy density in the reference frame of the emitting plasma 
by a factor of C 2 , and the mean photon energy by a factor of r. These two effects make it 
possible to explain the observed X-ray emission with electrons with modest electron Lorentz 
factors on the order of y = 100. The second model requires much higher Lorentz factors: 
assuming equipartition magnetic fields, electrons with Lorentz factors y ~ 10 7 are needed to 
explain the X-ray emission. Both models have severe shortcomings. For example, in the first 
model the long cooling times of the low-y electrons are at odds with the well-defined knotty 
structure of some X-ray bright jets. A weakness of the second model is that there is not yet 
a good explanation for the existence of two distinct electron populations with very different 
spectral properties. For lower-power FR-I-type sources like M87, Cen A, or 3C 31, the com- 
bined radio, optical and X-ray energy spectra are consistent with a synchrotron-origin of the 
X-ray emission from a single population of electrons. For an equipartition magnetic fields 
of B ~ 100 pG, the X-ray emitting electrons have Lorentz factors y ~ 10 7 , and radiative 
synchrotron cooling times of a few years. The X-ray emitting electrons thus do not move far 
from their acceleration sites before they loose their energy, and the X-ray brightness profiles 
tracks the acceleration of the high-energy electrons. The diffuse appearance of some jets 
implies quasi-continuous acceleration. 

A few radio galaxy jets have been studied with the Hubble Space Telescope giving not 
only high-resolution images of the optical brightness but also of the optical polarization. 
The magnetic field probed by the optically emitting electrons seems to be aligned parallel 
to the jet flow for the most part. However, upstream of the brightness maxima, the field are 
perpendicular to the jet flow. The radio and optical polarization behavior differs, indicating 
that the emission at different wavelengths samples different regions of the jet (Perlman et al. 
1999). 
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5 Acceleration to Ultra-high-Energies 

5 . 1 Some Properties of UHECRs 

One of the 11 fundamental science questions for the 21st century listed in the final report 
of the 2002 Decadal Review (Turner et al. 2002) is the nature of cosmic rays. The detection 
of cosmic rays at ultra-high energies (UHECR) dates back to the early Sixties, but only 
during the last 20 years detectors of sufficiently large size have become operational, that the 
origin of UHECRs can be addressed (for a review see Kotera and Olinto 2011). UHECRs 
are indirectly detected by observing the airshowers they trigger in the atmosphere. One 
can look for fluorescence emission or other radiation produced high in the atmosphere, or, 
alternatively, one registers the passage of secondary particle in charged-particle detectors 
on the ground. The AUGER observatory (Abraham et al. 2004) and the Telescope Array 
(Kawai et al. 2008) combine both techniques. 

The main observables used to infer the properties of UHECRs are the anisotropy, the 
composition, and the spectrum. The composition is difficult to determine, because the ap- 
pearance of giant airshowers can only be modeled with particle-physics event generators 
that involve extrapolations of behaviour observed in accelerator experiments, for the CoM 
energy of an arbitrary nucleus of 10 EeV energy with a nitrogen nucleus at rest is far higher 
than that achievable with even the largest man-made accelerator, the LHC at CERN. Con- 
siderable systematic uncertainty thus overshadows attempts to study the composition. It 
appears that between 1 PeV and 0.1 EeV we observe a trend from a predominance of light 
particles to heavy nuclei. Around 1 EeV, the composition is light again (Abbasi et al. 2010; 
Abraham et al. 2010b). Above 1 EeV, Auger observes a transition to heavier particles that is 
not see with other experiments at this time, possibly on account of statistics. 

The anisotropy is low around 1 EeV, where upper limits near 1 % have been published 
for the sidereal dipole anisotropy (Abreu et al. 2011). At higher energies above 57 EeV, for 
which little deflection would be expected, if the primary particles were protons, a correlation 
is observed between the arrival direction of particles and certain types of nearby AGN which 
in the end are proxies of the matter distribution within ~75 Mpc from us (Abraham et 
al. 2008; Abreu et al. 2010). The distance limitation is expected because nuclei at these 
energies undergo photodisintegration and photomeson production that provide losses on 
corresponding time scales, leading to the so-called GZK cut off. 

Whereas cosmic rays approximately obey a power-law spectrum with index s — 2.7 
( dN/dE oc E~ s ) below the so-called knee in the spectrum at 3 PeV, the spectrum of UHE- 
CRs is soft between a few PeV and 3 EeV with a power-law index s — 3 (Apel et al. 2009). 
At 3 EeV the spectrum hardens to s 22 2.6, a feature known as the ankle. Above 30 EeV one 
observes a flux suppression that has been identified with the GZK cut off (Abraham et al. 
2010a) A recent proposal (Aloisio et al. 2011) interprets those features as a proton cut off 
around 10 18 eV and another one around 3 x 10 19 eV associated with iron nuclei. 1 

5.2 Implications of the Maximum Energy 

Cosmic rays at energies below 1 PeV are almost certainly galactic in origin, and those at 
energies above 10 EeV are most likely extragalactic, but considerable uncertainty exists at 
intermediate energies. It is unclear at what energy the local cosmic rays turn from being 


1 This scenario implicitly postulates a very large proton-to-helium ratio in the source. 
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predominantly galactic to being mostly extragalactic. The relevance of this uncertainty for 
modeling the sources of cosmic rays is obvious: if the particles in the energy band above the 
knee at a few PeV, or above the iron knee at ~10 17 eV, are extragalactic, then considerable 
linetuning is required in matching the galactic and extragalactic components, because the 
spectrum softens at the knee. On the other hand, if cosmic rays up to a few EeV, i.e. up 
to the ankle in the spectrum, are galactic, then no such linetuning is required, but we need 
to identify the sources of EeV-band cosmic rays with objects present in the Galaxy, e.g. 
supernova remnants (SNR), pulsars, etc. This can be difficult, not because the source in 
question would not accelerate particles to high energies, which in fact we observe happening 
in SNR and pulsars, but because it is questionable that EeV energies can be reached. 

In fact, for typical interstellar magnetic field values, SNR shock fronts can hardly accel- 
erate cosmic rays to a PeV (Lagage and Cesarsky 1983a, 1983b). Particle confinement near 
the shock is supported by self-generated magnetic turbulence ahead of and behind the shock. 
Various plasma instabilities driven by cosmic rays can contribute to excite the turbulence to 
high levels, although which dominates remains an active topic of research. In the case of 
SNRs, that which has received most attention so far is the so-called streaming instability 
seeded by the cosmic -ray net current (e.g. Wentzel 1974; Skilling 1975; Achterberg 1983; 
Lucek and Bell 2000; Bell and Lucek 2001), and more recently its non-resonant counter- 
part (Bell 2004, 2005; Pelletier et al. 2006). In contrast, relativistic shocks operating at the 
interface between AGN/GRB flows and the surrounding medium reveal a short precursor, 
which restricts the plasma instabilities to small scale modes (Medvedev and Loeb 1999; 
Pohl and Schlickeiser 2000; Pohl et al. 2002; Reville et al. 2006; Pelletier et al. 2009; 
Lemoine and Pelletier 2010, 2011a), as discussed in detail in the following. Clearly, 
the amplitude of the turbulence sets up the pitch-angle scattering frequency and thus 
the acceleration rate (Maikov and Diamond 2001). In addition, it also sets the scale 
for the maximum energy, to which a remnant may accelerate particles. Although ana- 
lytical and numerical estimates suggest that cosmic rays can very efficiently drive mag- 
netic turbulence ahead of the shock (e.g. Bykov et al. 2011), so the turbulent magnetic 
field may be much larger than the homogeneous interstellar field (Lucek and Bell 2000; 
Bell and Lucek 2001), large increases in the magnetic field strength do not necessarily trans- 
late into a significant increase in the maximum particle energy (Vladimirov et al. 2006). 

For relativistic sources such as AGN or GRB, relativistic shock acceleration can 
be invoked, but even there certain limitations arise (e.g. Gallant and Achterberg 1999; 
Achterberg et al. 2001; Pelletier et al. 2009; Lemoine and Pelletier 201 lb; Eichler and Pohl 
2011). To be shock accelerated, a particle that has crossed the shock toward the upstream 
must be overtaken again by the shock. Assuming the shock moves at Lorentz factor r s , the 
particle must have been deflected (by gyration or scattering) through an angle A 9 > 1 / -Tv 
while residing upstream. This deflection must be accomplished within a time At at least of 
order R$/ fisc, where R s is the shock radius at which the shock once again overtakes the 
particle. The factor 1 / f} s in A 0 arises for subrelativistic shocks on account of the small in- 
cremental energy gain per shock crossing. The particle must cross the shock ~ 1 //3 S times to 
double its energy, with each crossing requiring at least a significant fraction of a gyroperiod. 

The fastest possible deflection is provided by undisturbed gyration in magnetic field ori- 
ented perpendicular to the shock normal, for which the angle between particle momentum 
and shock normal, 0, increases linearly with time, and the deflection rate A0/ At must then 
obey A0/ At = /3c/r g > c/(R s r s ). This inequality sets a maximum energy £ max to which 
a particle can be accelerated, because the rate of change of angle presumably decreases with 
particle energy. As r g = pc/ZeB = fiE/ZeB, this corresponds to a maximum energy, in 
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the limit of relativistic particles, of 


E mm = ZeBR s r s (9) 

Such regular deflection only occurs in a magnetic field that is coherent over scales larger than 
the path length of the particle. In the short precursor of a relativistic shock, this restricts B to 
the background, undisturbed magnetic- field value, and it thus limits the maximal energy to a 
rather small value for typical interstellar-medium conditions. To be noted is that Eq. (9) may 
provide a rather academic limit, because in a relativistic shock a perpendicular magnetic field 
in the downstream region renders acceleration very inefficient. Also note that the magnetic- 
field strength, B, is supposed to be measured in the upstream frame of the shock, i.e. it is 
the ambient field in the source frame in, e.g., an external shock of a GRB. 

Note that scatter- free gyration cannot in general confine a CR particle to a subrelativistic 
blast wave in all three dimensions. The particle generally drifts off to the side after gaining 
the potential difference /3 S EBR S in energy. Some scattering is required which will reduce 
E max ■ Also note that we have neglected both adiabatic losses and drift to the periphery of 
the shock front and assumed that being overtaken by a spherical blast wave is sufficient for 
further acceleration. 

For relativistic shocks, escape through the lateral boundaries does not provide a stringent 
constraint on the maximal acceleration energy unless sideways expansion of the blast takes 
place: as viewed in the shock front rest frame, the particle is confined if its gyration radius 
r g,o|sh < B , with R_ the lateral extension of the shock front. Since r g o|sh — r g ,o /E g , with 
r g the upstream gyroradius in the background field, confinement leads to E max < Ty R± e B , 
which is not as restrictive as the previous expression if R± > Rs/E s . 1° that limit, sideways 
expansion of the blast is negligible and the overall dynamics resembles that of a spherical 
blast wave. 

What of often invoked magnetic-field amplification by cosmic-ray induced instabilities? 
The growth of plasma instabilities in the precursors of shocks is inevitable, and therefore in 
a realistic situation we cannot expect undisturbed gyration in perpendicular magnetic-field. 
In fact, random scattering is required if the large-scale magnetic field is oriented parallel 
to the shock, because otherwise particles could not return to the shock. Therefore, random 
scattering in small-scale fields will make acceleration at parallel shocks faster, and thereby 
increase the maximum energy, in particular at nonrelativistic shocks. 

The scattering mean free path can be written as A. ~ c/ Dee , where the angular diffu- 
sion coefficient is given by D e e = SO 2 /St ~ (e B rms / Tine) 2 1 / fie where St ~ / / fie is the 
scattering coherence time, over which the particle scatters by an angle SO, and 1 is the co- 
herence length of the magnetic field (e.g. Eichler and Pohl 2011; Plotnikov et al. 2011). 
At sub-relativistic shock waves, one must now impose r 2 /l < PsRs, with r g = pc/ZeB rms . 
For relativistic shocks, the condition for the particle to suffer a rms deflection l/T s over a 
timescale Rs/c reads r 2 /l < T^R S , so that these two equations can be combined into 

Emm — ZeB^slRsETs ( 10 ) 

which is less than the previous expression when / <^ir g . Such small-scale fields are expected 
in the precursor of relativistic shocks, which cannot exceed r g0 /r 2 , although there is then 
an ambiguity related to the reference frame of the small-scale magnetic inhomogeneities; for 
simplicity, we have assumed here that these magnetic inhomogeneities are at rest in the up- 
stream plasma. Thus, provided some large-scale perpendicular magnetic field exists, simply 
tangling the field on small scales, /, does not necessarily raise E max . Magnetic-field ampli- 
fication enhances the maximum energy only if it increases B 2 ms l. Note that the expression 
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E max = ZeBR s , often taken from Fig. 1 of Frill as (1984), is consistent with Eqs. (9) or (10) 
only if /3 S and r s are both of order unity. We stress again that the limits described here 
may not be reached at a real shock. Leakage from the precursor and the conditions down- 
stream must also be considered when evaluating the maximum energy and the acceleration 
efficiency. 

Equations (9) and (10) suggest that the sources of UHECRs are likely systems involving 
relativistic shocks. 

5.3 Sources of UHECRs 

Besides reaching the required particle energy, the sources of UHECR must also be powerful 
enough to provide the source luminosity needed to sustain the local flux of UHECR. Possible 
source candidates of UHECRs are active galactic nuclei (AGNs) (Biermann and Strittmatter 
1987; Takahara 1990; Rachen and Biermann 1993; Pe’Er et al. 2009), clusters of galaxies 
(Kang et al. 1997), Magnetars (Arons 2003; Murase et al. 2009), and gamma-ray bursts 
(GRBs) (Waxman 1995; Vietri 1995; Murase et al. 2006). Depending on the model, these 
sources may also dominate the energy range around 10 18 eV (Berezinsky et al. 2006; Wang 
et al. 2007; Murase et al. 2008). 

One can constrain the acceleration capabilities of various sources through the magnetic 
luminosity of these sources (e.g. Norman et al. 1995; Waxman 2005; Lemoine and Waxman 
2009), as follows. Let us assume that acceleration takes place in an outflow at radius r mov- 
ing with possibly relativistic velocity fi (and Lorentz factor C) towards the observer, so as 
to benefit from Lorentz boosting. We assume that acceleration proceeds with an accelera- 
tion timescale f acc = Ar g /c in the comoving frame, with A > 1. Then the maximal energy 
at acceleration is at least bounded by the condition t acc < r /(T fie), which means that the 
acceleration timescale must be shorter than the comoving age of the outflow. This limit can 
be rewritten in terms of the maximal energy in the observer frame, E m ax and in terms of the 
magnetic luminosity of the source, L B = r 2 0 2 r 2 ficB 2 /4 as calculated in the source rest 
frame in terms of the jet half opening angle 0 and comoving magnetic-field strength B \ 

£max< 10 20 eV A- x r- x 0- x r y2 ZL x £ 5 (11) 

with L , £,45 = L b / 10 45 erg/s. One can check that this bound remains robust in the small 
0 limit, meaning 0T — > 0 for which side escape becomes important, and in the small 
fi limit. This bound indicates that rather extraordinary luminosities are required to accel- 
erate particles to ultra-high energies, under rather general conditions, although the bound 
depends on the charge of the particle. For instance, if one derives the magnetic luminosity 
of blazars through a leptonic modelling of the spectral energy distributions, one concludes 
that only the rare flat spectrum radio quasars with jet powers > 10 44 . . . 10 46 erg/s can ac- 
celerate protons to ~10 20 eV, while other B1 Lac and TeV blazars (FR I analogs) with jet 
powers ~10 40 . . . 10 44 erg/s appear limited to ~10 18 -10 19 eV (Lemoine and Waxman 2009). 
From this point of view, more compact sources such as GRBs and magnetars appear favored. 
For instance, a GRB of apparent isotropic luminosity 10 52 erg/s with U ~ 1 00 may produce 
particles with energy as high as Z x 10 21 eV for a magnetic conversion factor If B = 0.01. 
As discussed in Sect. 6.3, mildly or sub-relativistic shocks in a relativistic flow are more 
efficient accelerators of protons than ultra-relativistic shocks and are excellent candidates 
for being sources of UHECRs, owing to the magnetic-field amplification at shocks. 

The paucity of FR 2 radio-galaxies in the GZK sphere (radius ~100 Mpc) capable of 
accelerating protons to ultra-high energies might be compensated by the acceleration of 
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heavier nuclei in the less powerful and more numerous FR I radio-galaxies. In particular, 
Ptuskin and collaborators have shown that if radio-galaxies inject of light to heavy elements 
with a rigidity dependent maximal energy following Eq. (11), L B being related to the radio 
luminosity, accounting for the radio-luminosity function, one could explain rather satisfac- 
torily the observed spectrum (Ptuskin et al. 2011). It is also intriguing that the Pierre Auger 
Observatory reports an excess of events in the direction of the nearby radio-galaxy Cen A 
(although this latter happens to lie in front of one of the largest concentrations of matter in 
the GZK sphere, the Centaurus supercluster). However, it would be very difficult to under- 
stand the observed pattern of anisotropy if one assumes that the highest energy particles are 
heavier than hydrogen in such scenarios (Lemoine and Waxman 2009). 

Besides the actual source physics, the evolution of sources, the number of accelerators 
within a source (Aloisio et al. 2007), and the variation of source properties (Kachelriess and 
Semikoz 2006; Berezhko 2008) will also shape the local spectrum of UHECRs. 

Many properties of UHECRs can be impacted by their propagation in intergalactic space, 
such as their composition through photo-disintegration or their spectrum through cascading 
via photo-meson and photo-pair production, but it is difficult to disentangle the propagation 
effects from the results of physical processes operating inside the sources of these particles. 

Estimating the source luminosity using observed quantities and the known population 
statistics of the sources in question is subject to considerable uncertainties. As an example, 
Eichler et al. (2010) have recently estimated the local UHECR source luminosity, assuming 
all particles above the ankle at 4 EeV are extragalactic, and compared that with the observed 
gamma-ray production rate of all GRB. In contrast to earlier studies (Waxman 2004; Le and 
Dermer 2007), not the MeV-band gamma-ray fluence was used, which likely represents a 
thermal pool, but the GeV-band emission observed with Fermi-LAT, which measures the 
non-thermal tail of the energy distribution in the GRB primary charged particles, which, if 
hadronic, is the part that could contribute to the UHECR flux. It turns out that the UHECR 
source luminosity is more than a hundred times higher than the total GeV-band photon 
output, which places severe constraints on UHECR models involving GRB. If one posits 
that the Galactic to extra-Galactic transition takes place at ~ 10 19 eV and the MeV gamma- 
ray fluence traces the nonthermal particle population, the particle output of GRBs is more 
commensurate with their photon output (Waxman 2010). 

For each source class, one can also estimate a luminosity function, that is the differential 
source density needed to integrate the contribution of the sources over cosmological red- 
shift. While only nearby (< 200 Mpc) sources may actually contribute to locally observable 
GZK-scale UHECRs, the interaction products of the particles from all more distant sources 
will feed a cascade of energy that is eventually observable as a component of the extra- 
galactic gamma-ray and neutrino background, which are two other cosmic messengers that 
are complementary to the charged particles. The former has been recently measured with 
unprecedented sensitivity up to 100 GeV (Abdo et al. 2010), and thus provides invaluable 
constraints on, e.g., the so-called dip models, which assume that essentially all particle above 
about 1 EeV are protons. The redshift of the onset of photopair production with the CMB 
would then naturally lead to an ankle at the energy where it is indeed observed (Berezinsky 
et al. 2006). More precisely, it is the cosmic evolution of the source class in question that de- 
termines how much energy is fed into an electromagnetic cascade and eventually reappears 
in the GeV-band background radiation, relative to the UHECR energy flux at the ankle. 

5.4 The Transition from Galactic to Extragalactic Origin 

An open problem in cosmic-ray astrophysics is at what energy we observe the transition 
from a Galactic to an extragalactic origin of particles. The limit on inferred source power 
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per unit baryon mass required to sustain Galactic UHECR in the [4 — 40] EeV range that is 
imposed by the observed anisotropy limits is smaller by nearly 3 orders of magnitude than 
what is required for an extragalactic origin, as calculated in Eichler et al. (2010), and it cor- 
responds to the power per unit mass of gamma rays from GRB (Eichler and Pohl 2011). This 
is not only confirmation of the hypothesis that UHECR beyond the ankle are extragalactic, 
it also suggest that their sources are systems not persistently present in the Galaxy. Any as- 
trophysical source class, that is capable of accelerating particles to very high energies and 
should exist in the Galaxy, may fall short of accounting for the trans-ankle UHECR, but may 
nevertheless significantly contribute to the observed cosmic-ray flux between the knee and 
the ankle. For example, the numerical coincidence fits the hypothesis of a GRB origin for 
the Galactic component of UHECR, without invoking a much larger unseen energy reservoir 
for GRB. 

The interesting question is the rate with which such sources appear in normal galaxies 
such as the Milky Way. In other words, what is the role of intermittency? Generally, GRBs in 
the Galaxy are expected every million years or so, the exact rate depending on the beaming 
fraction and the detailed scaling of long GRB with star formation and metallicity. Therefore, 
only a small number of GRB can contribute to the particle flux at the solar circle, and their 
relative contribution depends on the location and explosion time of the GRB. Variations in 
the local particle flux must be expected, and neither the particle spectrum from an individual 
GRB nor the spectrum calculated for a homogeneous source distribution are good proxies. 

Pohl and Eichler (2011) have calculated the time-dependent transport of UHECR in the 
Galaxy, assuming it can be described as isotropic diffusion. They find that intermittency be- 
comes serious if the mean free path for scattering exceeds 100 pc, unless the source rate is 
much higher than 1 per Myr. On average, Galactic long GRB need to contribute only about 
10 37 erg/s in accelerated particles to fully account for the observed particle flux at 10 18 eV, 
assuming a Bohmian mean free path at this energy. UHECR from Galactic long GRB can 
meet the observational limits on anisotropy only if the mean free path for scattering is suffi- 
ciently small. Contributing the observed sub-ankle particles (at 10 18 eV) requires Bohmian 
diffusion if the UHECR are as heavy as carbon. A light composition such as protons or 
helium requires sub-Bohmian diffusion, which is a highly unlikely situation for isotropic 
diffusion. 

Much of the UHECR anisotropy arises from the expected location of long GRB in the 
inner Galaxy. Observations of GRB host galaxies suggest that regions of low metallicity 
and active star formation may be the preferred sites of long GRB (Levesque et al. 2010; 
Levesque 2011), which may skew the galactocentric distribution of long GRB toward the 
outer Galaxy. As there is no power problem with Galactic GRB, it may be worthwhile to 
also consider short GRB . They provide supposedly less power as a population, but they may 
have a very extended spatial distribution in the Galaxy (Berger 2010), leading to a reduced, 
but on account of intermittency not disappearing anisotropy. 

These conclusions can be applied with little change to the case of an origin of UHECRs 
in SNRs, assuming very efficient magnetic-field amplification can increase their ability to 
accelerate particles to energies significantly higher than 1 PeV (e.g. Ptuskin et al. 2010). The 
spatial distribution in the Galaxy of long GRB and SNR can be expected to be similar, and 
therefore the average anisotropy is the same for both long GRB and SNR. If one combines 
such a galactic component with a dip model, so that the galactic/extragalactic transition oc- 
curs below 1 EeV, Bohm diffusion and a mixed composition of the Galactic component may 
still be viable, given the systematic uncertainties in the measurements. It would be highly 
desirable to improve anisotropy measurement between 0.1 EeV and 1 EeV, and likewise 
better constrain the composition. 
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5.5 UHECR Summary 

Recent progress in UHECR research has built on data from new large-scale observatories. 
The interpretation of measurements of the composition, anisotropy, and spectrum of par- 
ticles provides constraining links between these observables, that are further strengthened 
by new precision measurement of, e.g., extragalactic gamma-ray background emission. The 
very low anisotropy observed for EeV-scale particles provides a strong limit on the con- 
tribution of Galactic sources, if the composition is indeed light as suggested by data. The 
anisotropy found above 60 EeV would be difficult to understand if the particles were heavy, 
which is suggested by Auger data, but not HiRes. If these particles were light, dip models 
might be favorable which, however, must be carefully constructed to not overproduce the 
50-GeV-scale gamma-ray background. 

The main obstacle to further progress clearly is the systematic uncertainty arising from 
the interpretation of the evolution of giant airshowers with particle-physics models that are 
extrapolated over at least 1.5 decades in CoM energy from the range testable with manmade 
accelerator experiments. 


6 Particle Acceleration at Relativistic Shocks 

Strong shocks occurring in astrophysical flows often generate power-law distributions of 
very-high-energy particles. This is the origin of most high-energy phenomena in astro- 
physics. The favored mechanism for the generation of supra-thermal particles is the famous 
Fermi process. It involves with the scattering of high-energy particles off magnetic distur- 
bances that allow them to cross the shock back and forth and thus to gain energy. Many 
studies have been performed in the 80 s and 90 s by assuming pre-existing magnetic turbu- 
lence. However, it turns out that the pre-existing turbulence is generally not strong enough 
to account for the acceleration performance. The nonthermal X-ray emission from SNRs 
(Cassam-Chenai et al. 2004), but see also Pohl et al. (2005), revealed that the magnetic field 
is strongly amplified in the vicinity of the forward shock. Recent theoretical studies have 
shown that the penetration of accelerated particles in the shock upstream flow can gener- 
ate magnetic turbulence that reaches a level much larger than the intensity of the ambient 
mean field (Bell 2004; Pelletier et al. 2006). In producing turbulence the cosmic rays loose 
a fraction of its global energy (about 10 percent of the incoming energy) but increases the 
maximum energy of particles (cf. Sect. 5.2). The turbulent field can reach an intensity of a 
few hundreds of pG, much larger than the value of a few pG of the ambient magnetic field 
in the Galaxy. 

These results incited similar investigations for relativistic shocks. Very encouraging re- 
sults were obtained around the turn of the century which extended the theory of Fermi 
process to the case of relativistic shocks and predicted the formation of a power-law en- 
ergy spectrum with an index .s = 2. 2-2. 3 and an acceleration time as fast as the Lar- 
mor time (Bednarz and Ostrowski 1998; Gallant and Achterberg 1999; Kirk et al. 2000; 
Achterberg et al. 2001; Ellison and Double 2002; Lemoine and Pelletier 2003). But disap- 
pointment came once the effect of the ambient magnetic field had been taken into account, 
because it inhibits the Fermi process even when one considers a strong Kolmogorov turbu- 
lence (Niemiec et al. 2006; Lemoine et al. 2006). 

In the following, it will be shown how the paradigm of the three interdependent aspects 
of collisionless-shock physics successfully works in the absence of any mean field: structure 
with a partial reflection on a barrier, supra-thermal-particle generation, magnetic-turbulence 
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generation. Then the scattering issue in the presence of a mean magnetic field will be ad- 
dressed and the requirement for circumventing the inhibition effect will be stated. Then an 
unusual fact in astrophysics will be emphasized, namely the necessity of considering some 
unavoidable micro-physics, that turns out to be crucial not only for the relativistic shock for- 
mation but also for making the Fermi process operative and producing high energy particles. 

6.1 Successful Fermi Process at Very Low Magnetization 


The most favored process for the generation of supra-thermal power law distributions is 
the Fermi process at shocks. Under astrophysical conditions the plasma flow that experi- 
ences a shock is supposed to carry a frozen- in turbulent magnetic field which allows particle 
scattering, and thereby permits particles to gain energy at each Fermi cycle, i.e. a cycle 
upstream-downstream-upstream or downstream-upstream-downstream. 

At a non-relativistic shock of speed f3 s = V s /c <<C 1 , the average gain per cycle is small, 
G = 1 + | '-y-fis (where r is the compression ratio, that reaches the value 4 when the shock 
is adiabatic and strong). Flowever this is compensated by a large number of shock cross- 
ings; indeed the escape probability (i.e. the probability for a particle to be entrained by the 
downstream flow and to not come back to the shock front) is low, P esc = Afi s /r\ the return 
probability P ret is thus large. A power-law distribution of energy is set up with an index that 
is a simple function of the compression ratio, in the non-relativistic case: 


s = 1 — 


In Fret 
In G 



( 12 ) 


Strong adiabatic shocks provide a particle spectrum with an universal index, s — 2, which 
is modified by losses, radiation losses for the electrons, expansion or escape for protons. 
Subsequent to escape, the spectrum is then steepened by the effect of diffusive propagation 
and escape of particles from the Galaxy. 

A sizable fraction of the incoming energy flux is converted into cosmic ray pressure: 

Per = lcr PuV; with £ CT ~ 0.1. (13) 


The successive Fermi cycles produce a precursor of supra-thermal particles (mostly pro- 
tons) of large extension (the diffusion length increases with the particle energy) and this 
penetration in the upstream medium (the ambient medium for an external shock) triggers 
MHD turbulence through two types of streaming instability, one is resonant and has been 
considered for many years (see for instance McKenzie and Voelk 1982), the other is non- 
resonant and has been considered more recently (Bell 2004; Pelletier et al. 2006), as briefly 
discussed earlier. That latter case is quite interesting, first because it is a simple and robust 
mechanism based on the supplementary Lorentz force associated with the plasma current 
that compensates the cosmic-ray current, second because it leads to a turbulent field of large 
intensity; indeed this latter can become much larger than the ambient magnetic field. The 
theory indicates that the fraction of incoming energy flux converted into magnetic energy 
can reach ~ /3 S , which is a few percent in SNRs, where one defines 

^=$BPuV s 2 . (14) 

A very important remark is that the efficiency of the Fermi process depends on the efficiency 
of the scattering process. By the way, the mechanism of Fermi acceleration is a simple 
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Table 2 Comparison non-relativistic shocks and relativistic shocks 

At non-relativistic shocks 

At relativistic shocks 

Weak escape probability 
Many cycles of weak energy gain 
Power law distribution s~ s with ~ 2 
Upstream distribution weakly anisotropic 
Partial reflection at shock front 
Generation of MHD turbulence upstream 
Acceleration time t aC c ~ W ft 2 

Significant escape probability 
Few cycles of large energy gain 
Power law distribution e~' s with j ~ 2.3 
Upstream distribution strongly anisotropic 
Partial reflection at shock front 
Generation of e.m. micro-turbulence upstream 
Acceleration time t acc ~ z s 


process, but the scattering, that controls the efficiency of the acceleration process, is the 
main issue. 

As for relativistic shocks, there are similarities and some differences with the non- 
relativistic ones, as summarized in Table 2. There are strong arguments that there is a sig- 
nificant generation of magnetic turbulence at the external shock of a GRB (Li and Waxman 
2006) and there is an obvious power-law distribution of ultra-relativistic electrons that syn- 
chrotron radiate, with an index compatible with the theory of the Fermi process at ultra- 
relativistic shocks (s = 2.2-23). The ambient magnetic field is very low and at first ap- 
proximation can be neglected. A remarkable work was published by Spitkovsky (2008) 
that fully validates the paradigm, combining three fundamental processes: the formation 
of a collisionless relativistic shock front with reflected particles, the generation of mag- 
netic turbulence and the generation of a power-law distribution through the Fermi process. 
This is a PIC (Particles In Cell) simulation of the development of a collisionless shock in 
a pair plasma (electrons and positrons) that runs with a Lorentz factor r s of a few tens 
(r s = (1 — P 2 )~ 1 - 2 ). The flow of reflected particles interacts with the flow of passing par- 
ticles leading to streaming-type instabilities, and the Weibel branch of instability describes 
the formation of intense small-scale magnetic filaments. The relevant scale of the physics 
is the inertial length (or skin depth) 8 = —. The spatial growth of the magnetic micro- 
turbulence produces a partial reflection of the incoming particles, which allows the forma- 
tion of a shock front, and self-consistently, the reflected particles generate the required level 
of micro-turbulence. Similarly as the non-relativistic case, conversion parameters £ cr , can 
be defined in the ultra-relativistic case: 

Per = ^p u r 2 c 2 , (15) 

^ = $ B p u r?c 2 . (16) 

And the simulations indicate that £ cr ~ 0. 1 and ~ 1-10 %, similarly to the non-relativistic 
case. The supra-thermal spectrum obtained in the simulation is close to the theoretical pre- 
diction with an index s — 2.4. Similar results were obtained later with PIC simulation in- 
volving a plasma of electrons and ions of (10 . . . 100)m e (Sironi and Spitkovsky 2009). 

6.2 Opening Phase Space with Finite Magnetization 

Many astrophysical shocks form in a plasma having a significant magnetization. The physics 
becomes more complex with a finite ambient mean field; it is controlled by the important 
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“magnetization” parameter a : 


B lf _ Bl sin 2 0 B 
47 zp u r 2 c 2 Anp u c 2 


(17) 


where Z? 0 is the field measured in the upstream flow frame (generally the ambient field), and 
B, f is the transverse component of the mean field measured in the shock frame. Like in 
non-relativistic shocks, the angle of the field lines with respect to the shock normal is very 
important. But whereas most non-relativistic shocks are in the so-called “sub-luminal” con- 
figuration, i.e. that the angle 0 B is not too close to 90° and thus particles can flow along the 
field lines, in ultra-relativistic shocks, it suffices that the field angle 0 B be larger than l/T s 
to prevent the return of particles to the upstream region. A generic ultra-relativistic shock is 
thus “supra-luminal”, and the magnetic field in the front frame can be considered as almost 
perpendicular, because its transverse component is amplified by a factor r s . This field ori- 
entation is a serious hindrance for the development of Fermi cycles. Neglecting for the time 
being any scattering process in a putative turbulence superimposed on the background field, 
the particle kinematics can be described as follows. A particle that enters the downstream 
flow of speed c/3 is dragged by the frozen in magnetic field and cannot easily come back 
upstream; it can be shown that it can come back just one time (Lemoine et al. 2006). Once 
upstream, it eventually comes back downstream, but in a subset of phase space that does not 
allow it to make a second cycle. Now, it might be thought that a strong turbulence could pro- 
vide efficient scattering allowing it to make several cycles. Flowever, the typical interstellar 
turbulent field with a large-scale coherence length behaves like an ordered magnetic field 
for such particles, because their penetration length upstream (f p = m p c 2 / T s eBo, measured 
in the co-moving upstream frame) is much shorter than the coherence length of turbulence 
(Lemoine et al. 2006). In self-generated small scale turbulence, scattering might be efficient 
enough to trigger Fermi acceleration, see below. 

The coherence length t c is formally defined as the range of the field correlation using the 
self-correlation function, C(r). For an isotropic turbulent state we can write (it can easily be 
properly modified in the case of anisotropic turbulence): 

f* OO 

tc= C(r)dr- (18) 

Jo 

which can be expressed as an integral over the turbulence spectrum, and one finds that for 
a spectrum proportional to k~P , the correlation length corresponds to large wavelengths for 
1 < /I < 2, as is the case of a Kolmogorov spectrum; for ()</!< 1 , the coherence length is 
in the shortest-wavelengths part of the spectrum. 

Moreover, the expected duration of the cycle would be much shorter than the eddy turn- 
over time of large-scale vortices. The requirements for efficient scattering off magnetic tur- 
bulence are quite challenging (Pelletier et al. 2009), for not only the intensity of the turbulent 
field must be much larger than the mean field, but also the coherence length must be shorter 
than a Larmor radius. When a scattering process develops, phase space is opened for op- 
erating a Fermi process if the scattering frequency is larger than the Larmor pulsation in 
the mean field. Short-scale turbulence leads to a scattering frequency v s oc s 2 , whereas the 
Larmor pulsation co L oc e; thus the range of particle energies for which the phase space is 
unlocked and Fermi process operative, is such that s < £ S catt = Ze(B 2 /B 0 )l c . 

At high magnetization (say a > 0.03) the shock is formed by generation of an intense 
coherent wave through a Synchrotron Maser Instability due to a resonance with the loop 
of reflected particles (Hoshino and Arons 1991; Gallant et al. 1992). The electromagnetic 
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wave propagating downstream is damped by synchrotron resonance and produces a thermal 
distribution. The wave propagating upstream carries away a fraction ~ 0. 1 a of the incoming 
energy in the case of an e + ~e~ -plasma; in a p + -e~- plasma, an electrostatic wake field is 
generated that heats the electrons up to equipartition while slowing down protons (Hoshino 
and Arons 1991; Sironi and Spitkovsky 2011). The formation of a power-law distribution 
(Hoshino 2008) has not been confirmed as far as we know. The final word has not been 
given on these issues of course, because the simulations have been conducted so far in ID 
(Hoshino and Arons 1991; Gallant et al. 1992) or 2D (Sironi and Spitkovsky 2011) over a 
limited amount of time. What happens in a more realistic 3D simulation, the dimensionality 
of which should allow more efficient cross-field transport, or on longer timescales, remain 
to be seen (see Jones et al. 1998). 

According to Sironi and Spitkovsky (2011), at lower magnetization, nothing happens 
except the thermalization of protons (T p ~ 0.2 r s m p c 2 ), until the magnetization reaches a 
very low critical value at which the Fermi process starts. As the magnetization decreases, 
indeed the precursor length scale increases, to the point where plasma microinstabilities trig- 
gered by the suprathermal particle population self-generate a small scale turbulence that can 
sustain the Fermi process (Lemoine & Pelletier 2010). Actually, one needs a very low mag- 
netic field to obtain an upstream penetration length of supra-thermal particles large enough 
for having a significant interaction of those particles with the incoming plasma and hav- 
ing a growth of micro-instabilities. The Fermi process works with the magnetic component 
of micro-turbulence at the inertial scale ~ 8 = c/a> P j. In principle it starts at even smaller 
scale, the inertial scale of electrons, however electrons are efficiently heated by the electric 
component of micro-turbulence and then the precursor becomes composed of electrons and 
protons of similar relativistic mass, like a pair plasma. This is a very interesting outcome that 
simplifies the physics which rapidly evolves towards conditions similar to those occurring 
in a pair plasma. Thus the PIC simulations of pair plasma are also valuable to understand 
the physics of shocks in electron-proton plasmas. Then a distribution function displaying a 
thermal part and a supra-thermal part with a power law is obtained. 

The transition towards the Fermi process is determined by the micro-instabilities that can 
grow when the upstream penetration of reflected particles is long enough. The fastest insta- 
bilities (Buneman instability, Oblique Two-Stream instability, see Bret et al. 2004) seem to 
essentially pre-heat the incoming electrons almost up to equipartition with protons. How- 
ever, more simulations are necessary to clarify this important point. The generation of mag- 
netic micro-turbulence by the Weibel instability, which is also studied in laboratory experi- 
ments, is thought to be the main ingredient to form collisionless shocks and to produce the 
Fermi process; however this is also under study by PIC simulations. The generation of mag- 
netic micro-turbulence occurs when the magnetization parameter falls below the following 
critical value (Lemoine and Pelletier 201 la), as confirmed by numerical simulations (Sironi 
and Spitkovsky 2011): 

<7<CTcrit=%- (19) 

F s 

Numerical simulations show that the level reached by that Weibel turbulence is such that 

= 1-10 %, which insures shock formation and Fermi process. Then there exists a large 
energy range for particle scattering when a § J . 

6.3 The Micro-physics Aspect of GRB Termination Shocks 

The main issue with Fermi processes based on the scattering off micro-turbulence is that 
the scattering frequency decreases as E~ 2 . The performance of Fermi processes at non- 
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relativistic shocks is determined by the scattering off large-scale, say Kolmogorov, turbu- 
lence which is fairly slow (much slower than the Larmor pulsation in the mean field) but 
decreases only as E~ 1 ^ . Thus, if we compare the Fermi process at relativistic shocks with 
the process at non-relativistic shock, this is like the hare and the tortoise: the scattering, and 
thus the acceleration rhythm, at relativistic shocks is very fast at low energy and decreases 
rapidly as energy increases, whereas, at non-relativistic shocks, it is slow at low energy but 
continues at higher energies with a moderate decline of its efficiency. 

6.3.1 Electron Acceleration and Radiation 


The external shock that drives the afterglow emission of GRBs may give rise to an effi- 
cient acceleration of electrons if the external medium is weakly magnetized, for the reasons 
discussed previously. If electrons thermalize with protons (as reasonably expected), their 
temperature is already very high at the beginning of the afterglow: T e ~ T p — 0.2 r s m p c 2 , 
which corresponds to a few tens of GeV. Intense short-scale magnetic turbulence develops 
because the interstellar magnetization parameter is very low, a ~ 10 -9 , whereas the critical 
value <7 cr it ~ 10 -6 , with r s ~ 300. 

What kind of radiation can be expected in such small-scale field, much more intense than 
the mean field? This depends on a so-called “wiggler” parameter a: 


a = 


s.1/2 


m, 
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m P c z m P 


(20) 


This parameter measures the capability of the magnetic force to deflect a relativistic electron 
of Lorentz factor y by an angle 1 /y (which is the reason why y does not appear in the def- 
inition). If a > 1, then the magnetic field produces a single deflection of the electron in the 
emission cone of half angle 1 /y, whereas if a < 1 the electron can undergo several wiggles 
in the emission cone. When a is large, the emission behaves like normal synchrotron radia- 
tion in a mean field, except that there is no polarization. When a is small, the emission is of 
“jitter” type (Medvedev 2000). Thus the emission caused by shocked and accelerated elec- 
trons at a relativistic shock is “synchrotron-like”, and the analysis of the emitted spectrum 
provides a diagnostic of the magnetic turbulence. 

It is quite remarkable that there exists an almost universal energy limit for the electron 
radiating in the intense small scale field (in agreement with Kirk and Reville 2010): 
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The corresponding maximum energy for the photons emitted in the quasi-homogeneous field 
is 
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( 22 ) 


where oif is the fine structure constant. The account for magnetic fluctuations of scales 
larger than the synchrotron emission formation length results in the photon spectra extended 
beyond the limit given by Eq. (22) (see Bykov et al. 2012). Thus a single synchrotron- like 
spectrum extending up to several GeV, even possibly a few tens, can be expected and is 
in fact compatible with observations. So the performance of relativistic shocks for electron 
acceleration and radiation appears very satisfactory. The conversion factor into radiation is 
trad ~ %B°Tnor s {y 2 ), and at the beginning of the afterglow £ rad ~ ~ 1-10 %. 
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6.3.2 Proton Acceleration Limited by the Fast Decay of Scattering 

Protons are expected to be accelerated at least as efficiently as electrons at ultra-relativistic 
shock waves. However, as mentioned previously, the ultra-relativistic Fermi process appears 
unable to push protons up to energies in excess of ~10 17 -10 18 eV, because the scattering 
time and thus the acceleration time increase with E 2 in the self-generated turbulent field, or 
scale with E, but then in the background unamplified field. 

Using one or the other, one does not find numbers significantly different from the limit 
associated to the mean field discussed in Sect. 5.2 (Eq. (9)): E max = Zr s eB 0 Rs — Z x (0.3 • 
10 7 GeV). Thus, although an energy of order 10 16 eV is achieved, which is something, the 
result is far from reaching the UHE-range. 

Precise performances of mildly or sub-relativistic shocks are not yet known and require 
more numerical simulations. However, some reasonable estimates are permitted by extrap- 
olating what we know about the two extremes: non-relativistic and ultra-relativistic shocks. 
The main guess is that we can expect a magnetic-field amplification at shocks with a con- 
version factor = 1-10%, occurring in MHD regime without severe limitation due to the 
super-luminal configuration, especially for oblique internal shocks (termination shocks in 
the hot spots of FR2 jets might be super- luminal). These assumptions can be applied to in- 
ternal shocks of AGN jets (in particular in Blazars jets), and to internal shocks of GRBs, as 
already discussed in Sect. 5.3. 

6.4 Conclusion and Prospect 

The triangular dependence of collisionless shock structure with a reflecting barrier for a 
part of incoming particles, with generation of supra-thermal particles and the generation of 
magnetic turbulence is a successful paradigm that applies to astrophysical shocks, both non- 
relativistic and relativistic. Numerical and theoretical works are making significant progress 
for both understanding the physics and providing quantitative results useful for astrophysical 
investigations. This includes not only the spectral index and cut off of the distribution of 
accelerated particles, but also the efficiency factors for the conversion into cosmic rays, 
magnetic turbulence and radiation. We have seen only the beginning of this line of study, 
which requires more PIC simulations and new types of hybrid codes involving relativistic 
MHD coupled with PIC codes for cosmic rays. 

The new results that have already been obtained are important. First, the strong am- 
plification of the magnetic field at SNRs received theoretical and numerical support; the 
astrophysical consequences are interesting, especially for our understanding of the Galactic 
contribution of the cosmic-ray spectrum. Secondly, current state-of-the-art PIC simulations 
indicate that the Fermi process does not operate at ultra-relativistic shocks with magnetiza- 
tion of order unity, which is supposed to be a frequent situation in high-energy astrophysics, 
as for instance in FR2 hot spots, in blazars, in pulsar wind nebulae. Such simulations need 
however to be extended both to higher dimensionality and to larger space-time domains 
before a definite conclusion can be reached. In particular, the issue of the non stationarity 
and/or corrugation of the shock front in relativistic regime should be investigated. Also the 
role of magnetic reconnections in the shock vicinity is a very important new topic whose in- 
vestigation is just starting. Whether and how acceleration proceeds in the mildly-relativistic 
regime also remains open for study. Thirdly the radiation processes that operate in most 
high-energy astrophysical sources involve relativistic electrons scattered in an intense short- 
scale magnetic turbulence; this leads to a renewed interest in the radiation physics with a 
view to use it as a diagnostic of the magnetic turbulent state. 
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The new trend in these topics is the important role imputed to micro-physics phenomena, 
which have a direct astrophysical impact. These developments incite interest in several other 
communities, including space-plasma physics, laser-plasma physics, astroparticle, and high- 
energy astrophysics. We live in exciting times. 
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